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Abstract. A recent first— principles approach to the non— linear rheology of 
dense colloidal suspensions is evaluated and compared to simulation results of 
sheared systems close to their glass transitions. The predicted scenario of a 
universal transition of the structural dynamics between yielding of glasses and 
non-Newtonian (shear-thinning) fluid flow appears well obeyed, and calculations 
within simplified models rationalize the data over variations in shear rate and 
viscosity of up to 3 decades. 



1. Introduction 

The rheological properties of soft materials, such as colloidal dispersions, presumably 
originate in a number of physical mechanisms, like shear induced phase transitions, 
direct potential and hydrodynamic interactions, advection of fluctuations, and shear 
banding or localization among others; see e. g. the collection of papers in [1]. 
At higher particle concentrations, the non-linear rheology depends on how steady 
shearing interferes with solidification during glass formation. Recently, we developed 
a theory for the non-linear rheology of dense colloidal suspensions aimed at this point 
[2]. It describes how the structural dynamics is fiuidized by advection of density 
fluctuations, while hydrodynamic interactions, non-linear flow profiles and ordering 
phenomena are neglected. Computer simulation studies can ensure that the latter 
processes are absent and thus provide crucial tests of the presented scenario. In 
this contribution, theoretical calculations will be compared to Brownian dynamics 
simulations of hard spheres by Strating [3] — without adjustable parameter in principle 
— , and to molecular dynamics simulations of a sheared binary Lennard-Jones mixture 
by Berthier and Barrat [4]. 

2. Theory 

2.1. General aspects 

A system of Brownian particles is studied in a prescribed steady shear solvent flow 
with constant velocity gradient and shear rate 7. The equation of motion for the 
temporal evolution of the many-particle distribution function is known [5], and has 
been solved for hard spherical particles at low densities [6]. This model constitutes 
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a first microscopic approach to real dense colloidal suspensions, and may serve as 
a model sheared glassy fluid [7]. It considers the "Brownian part" of the viscosity 
only, which, in Stokesian Dynamics simulations, Foss and Brady found to dominate 
compared to the hydrodynamic one for small shear rates [8] 

While the (approximate) approach developed in [2] gives general steady state 
quantities (like the shear distorted static structure factor) and their time-dependent 
fluctuations close to glassy arrest, we will concentrate on the thermodynamic shear 
stress (7(7) and the connected shear viscosity 77(7) = cr/7 + rjoo] here 7700 is the 
viscosity of the background solvent. The equations of motion exhibit a glass transition 
bifurcation, around which asymptotic expansions capture the transition from shear- 
thinning fluid flow to solid-like yielding. With the separation parameter e denoting 
the (relative) distance from the transition, and t a time scale obtained by matching 
onto microscopic transient motion, the following behaviors of a in the "structural 
window" have been established [2] 



where the Ci are positive material-dependent parameters (for the exponents 7 & a 
see e.g. [9]). The "structural window", here, is defined as the double regime |e| <C 1 
and |7io| -C 1, where the slowing-down of the structural dynamics dominates the 
steady state stress. While the divergence of the Newtonian viscosity 770 = toH~ 7 ci 
(first line of Eq. (1)) upon approaching the transition, applies to the linear-response 
regime of a fluid (e < 0), and is known from mode coupling theory (see the references 
in [9]), the novel predictions close to and above (e > 0) the transition describe the 
universal non-linear response of glasses to steady shearing with rate 7. Importantly, 
a "dynamic yield stress" cr + (e) = 17(7 — > 0+,e > 0) is obtained because a finite stress 
has to be overcome in order to force the glass to yield even for vanishingly small shear 
rate; a + is connected to the constants C2 and C4 in Eq. (1). While the yield stress 
varies strongly with distance to the transition deep in the glass, at fixed parameters 
close to the transition, the stress increases from er + with a power-law in 7, where 
the material-dependent exponent m lies around 0.15 in the models studied below. 
The given asymptotes are only the leading orders for e — > and -jt — > 0, while 
corrections can be obtained systematically [2], or are included in model calculations 
to be presented below. 

The dominance of the structural dynamics in determining 17(7*0, e) entails that 
all exponents or constants are functions of the equilibrium structure factor S q alone, 
except for the time scale to which matches to shorter non-structural dynamics. 
Thus, hydrodynamic interactions or inertial terms only influence the value of to, 
which ideally could be determined from an analysis of the intermediate scattering 
functions of the system [9] . This result arises because the small-shear rate rheology of 
glassy suspensions is dominated by steric hindrance (the "cage-effect") which is not 
qualitatively affected by the properties of the solvent around the particles. It is in 
agreement with the findings in Stokesian dynamics simulations [8] that shear thinning 
is dominated by a decrease of the Brownian part of the viscosity. The elimination of 
particle forces in favour of the quiescent-state structure factor S q is an approximation 
of unknown quality in the present situation, but in part motivated by the consideration 
of small shear rates. 




(1) 
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2.2. Models and simplifications 

The equations of motion, from which a in Eq. (1) follows uniquely for a given S q , have 
not been solved yet. Two approximate models were presented and discussed in [2] 
and shall be used in the following. While the schematic -model only incorporates 
the competition of two effects (divergent structural relaxation times with increasing 
e and loss of memory induced by shearing), the semi-microscopic ISHSM combines a 
semi-quantitative description of a quiescent hard sphere colloidal dispersion [9] with 
an isotropically-averaged shear advection of density fluctuations. Both models only 
depend on two parameters which map onto e and 7 introduced in Eq. (1), and thus 
can be viewed as minimal models for the described scenario. 

A problem when analysing data using both models arises from the ratio c 2 /ci 
in Eq. (1), which has a simple physical meaning. It gives the ratio of yield stress to 
transverse elastic constant (viz. shear modulus Goo) of the glass at the critical point, 
cijc\ = dia^/G^, where the numerical constant Ci = 1.0 (1) for the ISHSM (F^p). 
This ratio can be interpreted as a critical yield strain. Both models underestimate 
the effect of shearing leading to cr+/G^ = 0.33 (0.34) for the ISHSM (F$), 
while experiments give values around 0.05 indicating smaller strains are necessary 
for yielding [10]. While the schematic F-^-model is not meant to quantitatively 
capture such ratios, this error in the ISHSM presumably arises from the oversimplified 
handling of the shear-induced anisotropy of density fluctuations. The ISHSM treats 
all directions equivalent to the vorticity direction that is perpendicular to the flow 
plane. Perhaps unsurprisingly this underestimates the effects of shearing. We correct 
for this error in an ad-hoc fashion by rescaling the shear-rate 7 when considering 
77(7). For the F^-model this procedure is rigorously equivalent to an adjustment of 
the ratio a^/G^. 

3. Results and comparison with simulation data 

Before applying Eq. (1) to a solution of colloidal hard spheres at packing fraction <f>, 
at first the latter's critical value, <f> c , entering in e — C((j> c — 4>)/<fi c (with C = 1.5 
[9]) needs to be determined. This is done by testing whether the divergence of 
the quiescent viscosity (and corresponding structural relaxation time) for e — > 0— 
is observed. The inset of fig. 1 shows viscosities from experiments [11] and from 
Brownian dynamics simulations [3]. Also included arc self-diffusion coefficients from 
[12], which are predicted to vanish with D oc |e| 7 . Replotting the data with the 
calculated 7 = 2.62 [9], fits to the data above <j> > 0.50 give 4> c = 0.57 [12], 0.58 
& 0.60 for D, r\ (two outliers neglected) and simulations, respectively. Interestingly, 
the two experimental data sets provide rather close estimates for (f> c and indicate a 
strong coupling of diffusion and viscosity, Drj — ► 0. ^-q^D^ for e — > (neglecting the 
difference in <f> c ), with Dq the dilute single particle diffusion coefficient. The numerical 
factors is about half the predicted value [9] . We speculate that the discrepancy of the 
extrapolation of the simulation results arises in part because the data are not fully in 
the asymptotic regime. 

With the quantitative knowledge of e, only the matching time to is required to 
analyse the steady state viscosities in the structural window using the ISHSM. We 
chose to obtain it via the full fitting procedure which consists in matching by eye the 
numerical solutions to the non-Newtonian viscosity data. In this way, to is mainly 
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Figure 1. Steady state viscosities (symbols) from Brownian dynamics simulations 
[3] versus Peclet number Peo = -yd? /Do for packing fractions <j> as labeled. Fits by 
eye to the data for <j> > 0.50 with the ISHSM for separation parameters — e = 0.014, 
0.058, 0.097, 0.139 & 0.174 are given as solid lines and extrapolate to (fr c = 0.59. 
The matching time to = 0.019d 3 rj oo /fcBT is obtained and the theoretical 7 
is rescaled to 0.257 as discussed in the text. The inset shows a rectification 
plot with predetermined exponent 7 = 2.62 of viscosities from experiments [11] 
(circles) and simulations [3] (squares), alongside self diffusion constants (crosses) 
from [12], versus packing fraction. Linear fits to the data above <j> > 0.50 give 
(w/iloo) -1 /" 1 = 1.2s with (f> c = 0.58 (two outliers neglected), (ri/rj 00 )~ 1 ^"' = 1.6e 
with <t> c = 0.60, and (Dl/Dq) 1 ^ = 0.8 e with <j> c = 0.57 respectively. 



determined by the increase of the Newtonian viscosity, because i] ~ t a/(jto) holds 
and v/dto) becomes independent of to in the fluid for vanishing shear rate. The main 
panel of fig. 1 shows rj from the Brownian dynamics simulations as function of the 
dimensionless Peclet number Peo = 'yd 2 /Do, which measures the effect of shearing 
relative to the time a single particle diffuses its diameter d. The fits by eye using 
the ISHSM are included for packing fractions close to the transition, <f> > 0.50. 
From the fits, mainly from the divergence of 770 given in Eq. (1), the matching 
time to = 0.019d 3 77oo/fcBT is estimated, and inclusion of corrections to asymptotic 
behavior in the ISHSM-fits shifts the glass transition packing fraction closer to the 
other determinations; <f> c — 0.59 follows from the e used in fig. 1. Note that the 
solvent viscosity is included in the theoretical curves, ij — 77^ + t cr /(7^o)t- In 
the shear-thinning region, the viscosity diminuishes and approaches a behavior like 
r) ~ &t H with strong corrections, though, masking the power-law [2]. Because of 

\ The ISHSM calculations provide a[e, 'ytrjjPeo) for all values of e and Peo, while Eq. (1) captures 
the asymptotic behavior for e — ► 0, 7*0 — ► & Peo — > 0. Because we aim at describing the proximity 
of the glass transition, we match the parameters of Eq. (1) (</> c and to) there. Without matching the 
ISHSM gives <j> c = 0.52 & t = 0.025d 3 r; tx >/fe B T. 
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Figure 2. Symbols are shear stress (main panel) and viscosity (inset) data of a 
super-cooled Lennard-Jones binary mixture in reduced units taken from Ref . [4] ; 
from top to bottom, the temperatures are 0.15, 0.3, 0.4, 0.45, 0.5, 0.525, 0.555 & 
0.6 while T c ~ 0.435. The solid lines give fits by eye using the Fjg model from 
Ref. [2] with separation parameters: e = 0.050, 0.037, 0.021, 0.0, -0.027, -0.042, 
—0.054 & —0.083 (from top to bottom); units are converted by a = 'yr] = 1.8(t}7, 
where 7 = 0.53 Pe , and (r) = (-r(e,Peo)) [2]. 



the overestimate of cr+ in the ISHSM, this decrease would set in at too high 7 values 
only. In order to correct for the quantitative error, the theoretical curves are plotted 
versus rescaled shear rate, 7 * 0.25; i.e. a/(jt ) = f v (j * 0.25). With this ad-hoc 
correction, satisfactory agreement of theory and simulation results is seen for Peo < 1, 
where the steady state viscosity varies over two orders on variation of shear rate and 
packing fraction. For larger Peclet numbers, the data presumably lies outside the 
structural window where Eq. (1) applies. Motivated by numerical findings in [2], 
we speculate that the enhanced, ^-dependent steady state viscosities around Peo = 
1 - 10 in the Brownian dynamics simulations originate in the hard core repulsion. 
If so, hydrodynamic interactions which prevent particles from close contact could 
appreciably affect ?/ in this region. 

A second set of steady state shear stresses and viscosities is provided by recent 
large scale molecular dynamics simulations of a sheared simple liquid (a binary 
Lennard-Jones mixture) above and below its glass transition temperature [4] . Because 
kinetic parameters do not enter the theoretical predictions, and as linear flow profiles 
were obtained in the simulations, the universal predictions of our approach can be 
compared again. Figure 2 shows stationary shear stresses from the simulations and 
fits by eye using the F^p -model as specified in [2]. The model provides a relaxation 
time (r) as function of e and a dimensionless shear rate, denoted as Peclet number 
Peo, which are mapped onto the data as specified in the figure caption. The data 
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nicely span the glass transition temperature, T c w 0.435 already known [4], and are 
well compatible with a transition from a shear-thinning fluid to a yielding glass with 
finite yield stresses at and below the transition temperature. 

4. Conclusions and outlook 

We presented results of a microscopic theory of the nonlinear rheology of colloidal 
fluids and glasses under steady shear [2], and compared them with simulation and 
experimental data. This brought out the existence of a universal transition between 
shear-thinning fluid flow, with diverging viscosity upon increasing the interactions, 
and solid yielding, with a yield stress that is finite at (and beyond) the glass point. 
Numerical calculations could explain simulation results over up to 3 decades in 
variation in shear rate and viscosity. A quantitative analysis of larger data sets is 
required in order to determine the theoretical parameters for both simulations more 
accurately than the estimates found here. 

The approach we outlined should be improved with respect to the handling of 
shear-induced anisotropics, and stress-induced effects. The latter may lead to shear 
thickening behaviour that, for many colloidal materials, occurs at higher flow rates 
than those addressed here. This avenue will be explored in a future paper [13] on a 
version of the schematic model which is modified to include explicit stress- (as well as 
strain-rate-) dependence. 
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